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I. INTRODUCTION 

Quantum information processing is the study of meth- 
ods and efRciency in storage, manipulation and conver- 
sion of information represented by quantum states. Many 
quantum information theoretic concepts are closely re- 
lated to geometry. For instance, trace distance and fi- 
delity, which come out of the study of distinguishability 
between quantum states, are closely linked with Bures 
and Fubini-Study metrics. (See, for example, Ref. |H for 
an overview.) A few quantum codes can be constructed 
by algebraic-geometric means And finding the opti- 
mal quantum circuit can be regarded as the problem of 
finding the shortest path between two points in a certain 
curved geometry 

Recently, a few metrics on unitary operators with 
quantum information applications were found. For ex- 
ample, Johnston and Kribs introduced the kth operator 
norm of an operator acting on a bipartite system by con- 
sidering the action of the operator on bipartite states 
with Schmidt rank less than or equal to k. The fcth op- 
erator norm can be used to study bound entanglement 
of Werner states as well as to construct several new en- 
tanglements witnesses d, Q . Rastegin studied the par- 
titioned trace distance which shares similar properties 
with the standard trace distance Q. The partitioned 
trace distance shines new light on exponential indistin- 
guishability and hence can be used to investigate certain 
quantum cryptographic problems Interestingly, both 
the kth operator norm and the partitioned trace distance 
are related to the Ky Fan norm 0, Q . 

By asking the question about the minimum resources 
needed to perform a unitary transformation (a question 
of quantum information processing favor), I report fam- 
ilies of related metrics and pseudo-metrics on the set of 
unitary matrices (a result of geometric nature) in this 
paper. In this regard, these matrices and pseudo-metrics 
are very different from trace distance, partitioned trace 
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distance and fcth operator norm as the latter are more 
closely related to a different quantum information pro- 
cessing problem, namely, the distinguishability of states 
and operators. 

An interesting consequence of the discovery of this 
new family of metrics on the set of unitary matrices 
is that it gives refined measures of the degree of non- 
commutativity between two unitary matrices beyond the 
standard yes-or-no answer. Remarkably, while quanti- 
fying the level of closeness between two quantum states 
can be done by tools such as trace distance and fidelity 
(see, for example, Refs. 1, 8] for an introduction on this 
matter) , little has been done on the quantification of the 
degree of non-commutativity between two unitary oper- 
ators. 

Actually, the metrics and pseudo-metrics reported in 
this paper are constructed from certain linear combina- 
tions of the absolute values of the arguments of eigenvalue 
of a unitary matrix. To prove that these constructions 
are indeed metrics and pseudo-metrics, specific inequal- 
ities concerning the absolute values of the arguments of 
eigenvalue of the unitary matrices U, V and UV have to 
be established. The precise statements to be proven can 
be found in Definitions [2] and |4] as well as Theorem [2] 
below. Note that various authors had shown the valid- 
ity of a similar inequality, with the weighted sum of the 
absolute values of the arguments of eigenvalue for each 
unitary matrix being replaced by the largest argument 
of the eigenvalue for that matrix, provided that an ad- 
ditional condition constrainingthe arguments of eigen- 
values of U and V is satisfied |9l-[l3|. However, it is not 
clear how to modify their proofs to show the validity of 
the inequalities stated in Theorem [5J In fact, only a few 
results on the relationship between arguments of eigen- 
values of unitary matrices U, V and UV are known. Most 
mathematical works along this general direction concen- 
trate on the study of spectral variations of normal matri- 
ces on the complex plane (rather than arguments of the 
eigenvalues) as well as relations between the eigenvalues 
of Hermitian matrices Hi , H2 and Hi + H2 (known as 
the Weyl's problem). (See, for example, Refs. [Tsl - flU 
for comprehensive surveys.) Besides, techniques used to 
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tackle the Weyl's problem, such as the min-max principle 
stated in Sec. III.l of Ref. [l^, cannot be easily adapted 
to the case of unitary matrices. In this regard, the proof 
of Theorem O is also of mathematical interest. 

I begin by asking what is the minimum resources 
needed to perform a unitary transformation in Sec. |lll 
Motivated by the result of this quantum information the- 
oretic question, I define two closely related families of 
binary operations on the set of unitary matrices dp{-, •) 

and ctp {■,■), where /2 is a real- valued vector satisfying 
a technical condition, in Sec. IIIII Then, in Sec. IIVI I 
prove that these two families of binary operations are 
families of metrics and pseudo-metrics on the set of all 
unitary matrices, respectively. In Sec. [Vj I apply the 
metric introduced in Sec. IIIII to measure the degree of 
non-commutativity between two unitary matrices. Fi- 
nally, I give a summary and outlook in Sec. I VI I 



Recently, I showed that given a time-independent 
Hamiltonian H and a normalized initial state = 
J2j ctj\(,j{H)), the time r needed to evolve |(/>) to a state 
in its orthogonal subspace satisfies 



T > 



(1) 



for any x G M. Here A is a universal positive constant 
independent of H and \(p). Note that J2j \<^j?\^]{H)-x\ 
is minimized when x equals the median energy of the 
system M. Recall that the median energy M satisfies 



E 



a, 



> 
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II. MINIMUM RESOURCES NEEDED TO 
PERFORM A UNITARY TRANSFORMATION 

Two commonly studied problems in the field of quan- 
tum information processing are the maximum efficiency 
of a particular quantum information processing operation 
as well as the minimum resources needed to carry out 
such operation allowed by the known laws of nature . 
Since a unitary operator U, whose eigenvalues can be 
written in the form e*^^ 's with 9j £ (— 7r,7r], is the result 
of time evolution of an Hamiltonian H, one may naively 
use the values of \9j\^s as indicators of the minimum re- 
sources (in terms of the product of the average energy of 
the system and the evolution time needed) required to 
implement U. Nonetheless, this idea has to be polished 
as there is no physical meaning for the reference energy 
level and the overall phase of a unitary operator has no 
effect when applied to a density matrix. 

To refine the above idea, I use the following result. 

Definition 1. Let H be an n- dimensional Hermitian 
matrix. I follow the notation in Ref. fldij by denoting 
the eigenvalues and singular values of H arranged in 



and 
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j: \](H)<M 
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So, M need not be unique; and any M obeying the above 
two equations can minimize |Q!jp|Aj(_ff) — x\. From 
Eq. ©, I get 



r > 



n 



(3) 



where ^E{H, is the so-called average absolute devi- 
ation from the median of the energy of the system. Since 
Eq. ([3]) turns out to be the best possible lower bound on 
T given only &E{H, |0)), I interpreted &E{H, as an 
indicator of the maximum quantum information process- 
ing rate of a system [l^ . 

The above result can be used to prove the theorem 
below which quantifies the resources needed to perform a 
unitary transformation. The detailed proof can be found 
in Appendix 



descending order by Xi(H)'s and s^(H)'s respectively. Tiieorem 1. Let {aj)^^^ be a sequence of complex num 

.... . J . . .' . I . X — ^T^, I 12 - II lO^i lO^ 



where the index f runs from 1 to n. I denote the normal- 
ized eigenvector of H with eigenvalue Aj(_ff) by |^|(_ff)). 



bers obeying Yl^=i ["^iP = 1 '^^'^ l^^iP ^ lo^aP > ■■■ > 
|q;„P. And let U be a given n x n unitary matrix. Then 



R = min 

2:e[0,27r) 



Ht 



mm 

cxp{—iHt/h)- 



max 

|0)6C(ff,(a,)) 



A&E{H, |(/)))t 



(4) 



exists, where C{H,{aj)) is the set of normalized state 
kets in the form '^ilCp(j)(-^)) some permuta- 

tion P o/{l,2, ••• Moreover, the H and that 

attain the extremum in Eq. ^ can be chosen to have 
median system energy and XjiH)t/h e ( — TT, tt] for all j. 



In particular, there exists a Hamiltonian H such that 
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Theorem [T] can be understood physically as follows. 
Recall that 2iE{H, |0)) is an indicator of the maximum 
quantum information processing rate. Since U — e"'^*/'' 
up to an overall phase, one may implement C/, for in- 
stance, by evolving the system with a slow quantum in- 
formation processing rate H for a long time or alterna- 
tively by evolving the system with a fast quantum infor- 
mation processing rate H for a short time. In this re- 
spect, it is the product of ^E{H, and evolution time 
t that characterizes the resources needed to implement 
U. More precisely, the value of R defined by Eq. ^ 
in Theorem [T] can be regarded as an indicator of the 
least amount of resources needed to carry out the unitary 
transformation U on the worst possible |0) in the form 
S '^j\^p{j)i^)) some permutation P of {1, 2, ... , n}. 
The larger the value of R, the more the average absolute 
deviation from the median of the energy of the system 
times the evolution time is needed to carry out the trans- 
formation on the worst possible 



\d\j{u-') - \e\j{u), (6b) 

and 

ej{u) < \0\^{u). (6c) 

Definition 3. Let a G M and U = 

e'^3^'^-'|0^(C/))((/)^([/)| be a unitary matrix. 
Then, U"" is defined to be the unitary matrix 

E,.e-^'(^)|0j({/))(0^^(C/)|. 

Definition 4. Let U, V be two n-dimensional unitary 
matrices. Let fl = {pi,fi2, ■ ■ ■ , fJ^n) with Hi > > 
• ■ ■ > > 0. / define 

n 

i=i 

dp{U,V) = u{UV-^)p, (8) 



III. DEFINING FAMILIES OF QUANTUM 
INFORMATION THEORY INSPIRED METRICS 
AND PSEUDO-METRICS ON THE SET OF 
UNITARY MATRICES 

Based on the quantum information theoretic analysis 
in Sec. |lT] and the fact that R in Eqs. Q and ([5]) are 
essentially a weighted sum of singular values of the oper- 
ator i log U ~ xl for some a; S M, I propose a number of 
measures of the minimum amount of the average absolute 
deviation from the median of the energy of the system 
times the evolution time required to implement a unitary 
operation. I begin by introducing a few notations. 

Definition 2. Let U be an n-dimensional unitary ma- 
trix. Generalizing the convention adopted in Refi [lA], 
I denote the arguments (from now on, all arguments 
in this paper are in principal values unless otherwise 
stated) of the eigenvalues of U arranged in descending 
order by Oj[U) 's (where the index j runs from 1 to n). 

That is to say, U = Ej e**^^'^^ where 
Sj{U) G (— 7r,7r] and \(f)j{U)) is a normalized eigenvec- 
tor of U with eigenvalue e'^J . Similarly, I denote 
by \d\^{U) 's the absolute values of the arguments of the 
eigenvalue of U arranged in descending order. Occasion- 
ally, I need to refer to the arguments of the eigenvalues 
of U without any specific order using the notation 9^ 's. 
And in this case, I denote the corresponding normalized 
eigenvector by \<p'j)- 



For example, let U 



Then, efiU) = 0, 



1 

-i 

e\{U) = -7r/2, \e\\{U) = 7r/2 and \e\\{U) = 0. 

Note further that for any n-dimensional unitary matrix 



W ) = 



(6a) 



i.(;7)Y =minz.(e"C/)^, 



and 



dY{U,V)^mmdpie"U,V). 



(9) 



(10) 



I am going to prove in Sec.|TV]that d^(-, •) and rfj(-, •) 
are a metric and pseudo-metric, respectively. This justi- 
fies the use of the notations. On the other hand, 
and are not a norm or pseudo-norm. This is be- 

cause U -\-V may not be unitary so that iy{U + V)i2 and 
h'{U + are ill-defined. Hence, it does not make sense 
to talk about validity of the (additive) triangle inequality 
for and i^(-)J- Nevertheless, I am going to prove 

in Sec. llVl that and i^(-)J satisfy the multiplicative 
triangle inequality. 

Note further that v{U)fi is a weighted sum of the sin- 
gular values of a certain infinitesimal generator of the 
unitary operator U . And from Theorem [TJ I know that 
;y([/)Y is indeed a measure of the minimum average ab- 
solute deviation from the median energy of the system 
times the evolution time needed to perform the unitary 
transformation U in the sense that 



x£[0,2t!) Ht: cxp{-iHt/h)=e"U 



max^£'(iJ, \(f)))t, 



where the maximum is taken over all state kets 
fJ- ^ .,ru\, 1^, .|2 



in the 

form E;=i«.ICp(,)W) with |a,f = and P 

is a permutation of {1, 2, . . . , n}. In contrast, iy(U)p can 

be regarded as an "unoptimized" version of v{Uy? in the 
sense that 



Ht: CKp{-iHt/h) = U 



-.&E{H,\4>)) t, (12) 
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where the maximum is taken over the same set as in 
the maximum in Eq. pT|) . In fact, the global phase of 
L/, which has no physical meaning, affects the value of 
y{U)fi but not j/(J7)J. Nevertheless, readers will find 
out in Sec. |V] that y{-)p, is useful in studying the degree 
of non-commutativity between two unitary operators. 

Three important points are stated. First, one should 
investigate the properties of Eqs. ([7))-([TU| for a general jl 
in order to get a more complete picture on the minimum 
resources needed to evolve different kind of initial states 
|(/))'s by the unitary operator U . 

Second, an inequality giving a lower evolution time 
bound like the one stated in Eq. (|3]) is sometimes known 
as a quantum speed limit ^6,]. In addition to Eq. ([3|), 
the other two important quantum speed limits are the 
so-called time-energy uncertainty relation r > ■Kh/2/S.E 
where AS is the standard deviation of the energy of the 
system and the Margolus-Levitin bound r > 7rh/2{E — 
Eq) where E and Eq are the average e nergy a nd ground 
state energy of the system, respectively One may 

define similar measures of the minimum resources needed 
to implement a unitary operator U based on these two 
quantum speed limits. Unfortunately, it is not very fruit- 
ful to pursue in this direction because these measures are 
not metrics. (Note that the absence of metric structures 
here only refers to the measure of minimum resources 
required to perform a unitary operation based on these 
quantum speed limits. It does not mean that the two 
quantum speed limits have no geometrical meaning. In 
fact, the time-energy uncertainty relation itself is closely 
related to the Bures metric [20j.) 

Third, is related to the distance between two 

quantum states as follows. Suppose l^*!) and |4'2) are 
two n-leveled pure states. Recall that the Bures angle be- 
tween them is X = |(4'i|^2)|- Then, amongst all the uni- 
tary transformations U obeying U\'^i) — |5'2), the one 
that minimizes '^{U)'? satisfies Oi{U) — x, Oj^{U) — ~x 

and dj{U) = for 1 < j < n. (One quick way to see this 
is that the required U only need to perform rotation on 
the Hilbert space spanned by and |4'2). Then the 
problem can be handled by standard minimizing tech- 
niques.) 



IV. PROPERTIES OF !/(-)m, d^(-,-)> AND 



and 



Since p = (mi, M2, • • • , M™) = (Mi - Mi+i)/^ 

^nfj^^^ where 



j entries 



at'-'I -(1,1,--- ,1,0,0,--- ,0), 
I conclude that 



(13) 



n-l 



j^({/)J > y^ip, - ti,+iHu)^^y, +tiMu)l^,. (14b) 



Furthermore, the coefiicients in the R.H.S. of the above 
two equations are non-negative. It turns out that one 
can deduce many properties of i'{-)p from the properties 
of '^{■)pu] 's. 

Here I list a few basic properties oiv{-)p, J^(-)Jj d,p{-, •) 
and c?Y(-, •). The proofs are left to the readers as they 
are straightforward. (Note that all J7's and V's appear 
below are n-dimensional unitary matrices.) 

• and iy'(-)'p are both functions from U{n), the 
set of all n X n unitary matrices, to [0, '!^Y^^=i Mjl- 
In fact, is a surjection with v[U)p = TrX^j 
for all fl if and only if [/ = —I. 

• v{U)p = and v{U)J = if and only iiU = I and 
U = e^^I for some a; G M, respectively. 

• z^(e"[/)Y = v{U)J for all x G R. 

. u{U-^)p = y{U)p - u{VUV-^)p and v{U'^)J = 
v{U)J^v{VUV-^)J. 

• i^{U)ap = av{\J)p and viJJ)^^^ = a.v{U)1 for all 
a > 0. 

• v{U'')p < \bW{U)p and iy{U^)'? < \bW{U)^ for all 
6 e M. Moreover, the equalities hold if |6||6'|t(J7) < 



. z.(C/)V , = 2KC/)^, and ^^(t/)^.,,, = ^(^)^.,, 
whenever 2j + 1 < n. 

• If [/ (t) is a continuous one-parameter family of uni- 
tary matrices, then v{U {t))p and v{U (t)) J are con- 
tinuous. This result is more involved. By the the- 
orem that roots of a poly nomial vary continuously 
with its coefficients |2ll . [2^ . eigenvalues of U{t) 
must vary continuously on the unit circle. And as 
the absolute value of argument of a complex- valued 
and nowhere zero function is continuous, iy(U{t))p 
is continuous. The continuity of v{U{t))'^ then fol- 
lows the fact that the pointwise minimum of a col- 
lection of continuous functions is also a continuous 
function if it exists. 

• If H{t) is an n-dimensional Hamiltonian and U{t) 
is the unitary operator generated by H{t), then 



dv{U{t))p 



dt 



(15) 



'^{U)p = ^{fJ'j - ^J'j+l)l'{U)pu] + ^ini^{U)pM (14a) 



This shows the relation between rate of change of 
v{U{t))p and the singular values of its generator 
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H{t). Actually, the R.H.S. of Eo^ 
called generalized spectral norm [2J 
tonian -ff(O). 



15)) is the so- 
of the Hamil- 



• dp{-, ■) and d'? {■, •) are positive definite and positive 
semi-definite functions, respectively. And they are 
both symmetric functions. 

• dp{U, V) — n J2j Mi for fl if and only if V = ~U . 

• Let Ui be rii-dimensional unitary matrices, then 

n2v{Ui)p_[r,^]+niv{U2)0r.2] and i^(J7i(g)C/2)^„i„2] < 
n2viUi)Y,„ , + nii/(C/2)Y„ ,. (Note that I have 
abused the notation a little bit as clearly the 
lengths of the vectors pL^'^^\ /^["■i"^] in the 

above inequalities are different even though each of 
them has the same number of non-zero entries.) 

Interestingly, except for the second property (the nec- 
essary and sufficient condition for v{-)fi and v{-)'^ to 
be zero) and the last property (concerning the tensor 
product of unitary matrices), all the above properties 
require only /ij > for all j. In contrast, the condition 
A^i > > • ■ ■ > ^ is required for dfi{-, ■) and dfi{-, ■) 
to satisfy the triangle inequality. In fact, the triangle in- 
equality follows directly from Eq. and the theorem 
below. And this theorem alone is also of great interest 
as it bounds certain weighted sums of the arguments of 
eigenvalues for the product of two unitary matrices. 

Theorem 2. Let U, V be two n- dimensional unitary ma- 
trices and jl satisfying the requirements stated in Defini- 
tion^ Then 

W{U)p - ,^{V)p\ < ,^{UV)p < v{U)p + v{V)p. (16) 



I move both the proof of Theorem [2] and the condi- 
tions for the second inequality in Eq. to become an 
equality to Appendix [B] as they are rather involved and 
technical. 

Corollary 1. Eq. \16\) in Theorem\^ also holds if v{-)p 
is replaced by i^{-)'p ■ 

Proof. For any n-dimensional unitary matrices U and 
V, there exist x,y e R such that v{e''''U)p = v{U)J 
and v{e^yV)p = v{V)J . By Theorem H v{U)J + 
v{V)Y > v{e'^='+y^UV)p > v{UV)Y. The proof of 



'^{V)'p < i^iUV)'? is then a replica of that 



of the first hah of Eq. (H 



□ 



From Theorem [5] and Corollary [1] I conclude that 
dp{-,-) and dY(-,-) are a metric and pseudo-metric, re- 
spectively. Furthermore, i/(-)^, i'(-)Y behave somewhat 



like a norm and semi-norm respectively in the sense that 
v{UV)p < i'{U)p -\- v{V)p. In other words, the addition 
operation in the conventional definition of a norm is re- 
placed by multiplication. Note further that v{ )p does 
not obey i'{aU)p — \a\i^{U)p for all scalars a. 

Remark 1. For any n-dimensional matrices U and V 
satisfying 9\{U) + e\{V) < n and OHU) + ei{V) > -n, 
the following closely related bound 



9i{uv) < ei{u) + 0i{v) 



(17) 



has bee n p roven by a few authors using different meth- 
ods 0, [71 . [71 /. Recently, our group J2A . [H/ also found 
elementary proofs of Eq. |_?7| ) as well as the special case 
of Theorem\^in which jl = jl^'^^ . 

Remark 2. The proof for the special case of fl = /^["l 
is rather simple. I simply need to use the identity 
YY^^^e'0]iuv)/2n ^ det[([/T/)i/2"] = [det(t/l/)]i/2" = 

[detJ/dett/]^/^" = det det 
n;^ie^K(^)+»|(^)]/2". Then, by equating the ar- 
guments on both sides and by observing that both 
arguments are in the interval (— 7r,7r], / obtain the 
required inequality. 

Remark 3. It is easy to check from the proof of Corol- 
laryUlthat iy{UV)Y = v{U)J +v{V)J if and only if there 

exist x,y e K such that v[u)J = i/(e"C/)/i, i'{V)J = 
v{e^W)p, v{e^'^='+y'^UV)p = i^(e"C/)/i + v{e'yV)p and 
iy{UV)Y = i^(e^''''+y^UV)p. Therefore, the conditions for 
equality of the triangle inequality for v{-)^ seems to be 
more stringent than those for v{-)p. On the other hand, 
as I have pointed out in the list of the basic properties 
ofv{-)J that KC/)J[3, = KC/)J[., = MU)l^- Thus, 
the conditions for v{UV)J^^^ = '^(t^)^H + '^(^)J[H 
not particularly more stringent than those for i'{-)p[m] for 
TO < 3. But in any case, as m or n increases, it is more 
and more difficult for v^UV)"? = J^(C^)Y + i'(^) J pi'o- 
vided that U and V are drawn randomly from the Haar 
measure of U{n). This finding is reflected in the plots 
in Fig. Interestingly, by comparing Figs. [7] and [H / 
find that on average ( over the Haar measure of U (n) ) 
the value of v{U)J -f v{V)J - v{UV)J is smaller than 
v(U)p -\- u{V)p — i>{UV)p. This is probably due to the 
fact that v{U)J < v(U)p. 



Remarkably, several new inequalities involving the 
eigenvalues of a product of two unitary matrices can be 
found by modifying the proof of Theorem [2] [1^ . But 
I will not further elaborate on this matter here as it is 
beyond the scope of this paper. 
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FIG. 1. Plots of i^{UV) against i^{U)^[m] + >-'{V)p[r„] for various values of rn where the dimension n of the unitary matrices 
U and V equals (a) n = 2, (b) n = 3 and (c) n = 4. The matrices U and V are randomly chosen from the Haar measure of U (n) 
and the sample size in each case is 1000. The black lines represent the situation in which t'(Z7y) -[m] = i'{U)p[m] + i^{V)p[m] ■ 
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V. MEASURE OF THE DEGREE OF 
NON-COMMUTATIVITY BETWEEN TWO 
UNITARY OPERATORS 

Although the family of j/(-)^'s are only "unoptimized" 
measures of the minimum resources needed to perform 
the unitary transformation in its argument, it is useful in 
measuring the degree of non-commutativity between two 
unitary operators. 

Definition 5. Let U, V two n- dimensional unitary ma- 
trices. Suppose fl be the vector defined in Def. / define 

GipiU, V) = dpiUV, VU) = iy{UVU-^V-^)p. (18) 



for i = 1,2. Then, 



d^€p{Ui{t),U2{t)) 



dt^ 



t=o 



21im 



€p{Ui{t),U2(t)) 

i2 



(19) 



This shows the relation between the curvature 
of <Lfi{Ui{t),U2{t)) and the singular values of the 
commutator of the the corresponding generators 
[H^{t),H2m. 



Physically, €p{U, V) can be interpreted as a measure of 
the minimum average absolute deviation from the median 
energy of the system times the evolution time required to 
convert UV to VU. The higher the required resources, 
the more "non-commutative" the operators U and V are. 

The following simple properties of £^(-,-) carry over 
from the properties of i^(-)/J reported in Sec. IIVI 

• Cp{U, V) : U{n) X U{n) — > [0, vr^J^^ ^ij] is a sur- 
jective map. 

• €p{U, V^) = if and only if U commutes with V. 

• Cp{-,-) is not a metric. In fact, any meaningful 
measure of commutativity is not a metric for com- 
mutativity is not transitive. 

• £^(e"C/,e*?^F) = Cp{U,V) for any x,y€R. 

>€piu-\v-^) = GipiU.V) 

<tp{wuw-^ ,wvw-^). 

• €ap{U, V) = ailp{U, V) for all a > 0. 

• €p{U,V) = TT^^./ij if UVU-^V-^ = In par- 
ticular, the three Pauli matrices a^, (Jy and CTz obey 
<^/7(o-x,o-j/) = ilp[ay,a^) = <tp{az,cr^) = Tr^^feMfe- 
In this sense, pairs of distinct Pauli matrices are 
examples of the most non-commutative unitary op- 
erators in U{2). 

• Vi(C/i®Vi, [72(^^2) < epi^{Ui,Vi)+€pw{U2,V2), 
and e:^[„i„,j {Ui ® Fl, C/2 Fa) < naC^-ji-ii (C^i, V^i) + 

^ie:^M(f/2,T^2). 

• If [/ (t) and V{t) arc continuous one-parameter fam- 
ilies of unitary matrices, then €p{U (t), V{t)) is con- 
tinuous. 

• Let Hi{t) be an n-dimensional Hamiltonian and 
Ui{t) be the unitary operator generated by Hi{t) 



VI. DISCUSSIONS 

In summary, I have introduced a family of metrics 
dp{-, •) and a family of pseudo-metrics dY •) on finite- 
dimensional unitary matrices. In particular, the pseudo- 
metric d'p {U,V) can be interpreted as a measure of the 
minimum average absolute deviation from the median en- 
ergy of the system times the evolution time needed to per- 
form the unitary transformation UV~^ (and hence equiv- 
alently also VU^^). Besides, dj(-, •) is related to the Bu- 
res angle between two quantum states; while dp{U, V) is 
related to the generalized spectral norm of the infinitesi- 
mal generator of UV~^ . 

Another aspect of this paper is the proposed measure 
on the degree of non-commutativity between two unitary 
matrices [/ and based on i'(f/F[/~^F~'^)^. Its physical 
meaning is a measure of the minimum resources needed to 
convert UV to VU . Interestingly, this measure is related 
to the generalized spectral norm of the commutator of 
the infinitesimal generators of U and V . 

The analysis here so far are restricted to finite- 
dimensional unitary matrices. It is instructive to see 
how it can be extended to cover the infinite-dimensional 
case. Also a possible future research direction is to in- 
vestigate the possibility of extending the results here to 
trace-preserving completely positive maps and to find 
non-trivial applications in quantum information process- 
ing. 

Appendix A: Proof of Theorem [T] 

The following lemma is needed to prove Theorem [T] 

Lemma 1. Let Xj 's and yj 's be real numbers satisfying 
xi > X2 > ■ ■ ■ > Xn and 2/i > 2/2 > • • • > J/n- Then 

n n 

"^Xjyp^j) <J2^3y3 (Al) 

i=i i=i 

for all permutation P of the set {1, 2, • • • , n}. 
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(a) 




01 23456789 10 11 12 
\\UV\\}» 

(b) 




1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 

wuvwy 

(C) 




2 4 6 8 10 12 14 16 



\\uv\\}» 

FIG. 2. Plots of v{UV)Y[^] against i^(?7)Y[„] + for various values of m where the dimension n of the unitary matrices 

U and V equals (a) n = 2, (b) n = 3 and (c) n = 4. The parameters and sampling method used are the same as those in Fig. [T] 
Since i^((7)Y[3] = = 2i^(r/)Y[i] 1 I only show the cases of m = 1 and m = 4 in this figure. And to allow easy comparison, 

the scales used in the plots are the same as those used in Fig. [T] 
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Proof. Eq. (jAip is trivially true for n = 1. And its valid- 
ity for n = 2 follows from the inequality (xi — X2){yi — 
2/2) > 0. This Lemma can then be proven by mathemat- 
ical induction on n. □ 

Proof of Theorem [H Note that t can only be when 
gtxjj _ J ^Yiis case, Eq. (U) clearly makes sense 
and is equal to 0. So, this theorem is trivially true. 



Note that the set {((^|/HTff |0) : |(/)) € C{H,{aj))} 
equals {Ej^i l^j ps^(j)(-ff) : -P G S',,}, where de- 
notes the permutation group of n elements. So from 
Lemma [TJ for any given H, t and x, the maximum in 
Eq. exists and is equal to ^iE"=i lajP^^^)/^- 

Now suppose s\{H) > ttS/^. I consider the new 
Hamiltonian H' formed by changing only the eigen- 
value A corresponding to s\{H) to A mod {27Th/t) = 
A' e {-Trh/t,Trh/t]. Clearly, |A'| < |A| = s\{H) and 



That is to say, v[iSLX\^j^ciH.{aj)) {<PW H\(j>) does not in- 
crease if H is replaced by H'. By repeating this proce- 
dure at most n times, I conclude that in order to look for 
the minimum in Eq. (jA2l) . I only need to consider those 
H's with all their eigenvalues in {—irh/t, Trh/t]. Note that 

the expression ma.x^^^^ciH,{o'j)){4'W H'^ H\4') depends on 
the eigenvalues of H rather than the eigenvectors of H . 
Moreover, for a fixed a; S K, there is a unique set of eigen- 
values Aj(7J)'s for H such that X\{H) G {-Trh/t, irh/t] 
and exp (—iHt/h) = e^^U . So, the second minimum ex- 
pression (that is, the minimum over Hi) in Eq. (jA2p ex- 
ists for each a; g K. As a; varies, the set of eigenvalues 
of H that minimizes the second minimum expression in 
Eq. (jA2l) changes linearly with x modulo 2Trh/t. Hence, 
this second minimum expression is a continuous function 
of a; e M with period 2tt. Therefore, R' exists. 

Now I go on to show that the H and that min- 
imize Eq. (|A2p can be chosen to have median system 
energy M = 0. First, I claim that l^jp > 1/2 where 



Whereas if i < 0, then 9E{H, \(p})t < 0. In this 
case, I may consider —H and —t instead in analyzing 
Eq. (HI) for &E{~H,\(j))){-t) > and exp{~iHt/h) = 
exp[~i{—H){—t)/h]. Hence, I only need to consider the 
remaining the case oi t > from now on. 

I first justify the use of maximum and minimum in 
Eq. ([?]); and I do so by considering a similar min-max 
expression: 



max m^m^ 

\(t>)eC(H,{aj)) h 

I 

ex-p{—iH't/h) — exp {~iHt/h). In other words, there 
exists fee {1,2, ... ,n} such that 

r sj^,{H) if J <fc, 

sj{H')^\ |A'| ifj^k, (A3) 
[sj{H) iij>k. 

Consequently, by Lemma [TJ 



(A4) 

I 

the sum is over ah j's with Xj{H)t/h e {tt} U (-7r,0]. 
Suppose the contrary, as H has at most n distinct eigen- 
values, there exits e < sufficiently close to such that 

{(t>\\J {H + e)^ {H + e)\(l)) < {<j)\\/ H\(l)) . This contra- 
dicts the assumption that H and | (j)) extremize Eq. (jA2l) . 
By a similar argument, ^ |ajp < 1/2 where the sum is 
over all j's with Xj{H)t/h £ [0, tt). Hence, fr om Eqs. (|2ap 
and (j2bl) , the median energy of the system M is equal to 
0. 

From the discussions between Eqs. ((T]) and ^ on 
the minimization of J2j I'^jl'^l'^ji-^') ~ I conclude 

that &E{H', 10)) < max|^)gc(ff',K))(</>|Vi/'ti/'|0) with 
equality holds if the median system energy M 0. 
By comparing Eq. ^ with Eq. (jA2p . I deduce that 
R< R' a R exists. Recall that Eq. (|X2|) is well-defined 
and the extremum can be attained by H and \(f)) such 
that Xj{H) g {—Trh/t, irh/t] and the median system en- 
ergy M = 0. With these H and &E{H,\(j))) = 



R = min min 

xe[0,2TT) Ht: Gxp(-i//t/fi)=e"C/ 



fc-1 



max {^'\Vh^\<I,') - max {ct>\VHTH\^) ^ \akf\X'\ - \a,\h\{H) + {\ai\ 

p')eC{H' ,{a,}) \<p)£C(H,(aj)) ^ 



< si{H) 

= 0. 



|afcP-|aiP 



fc-i 

Ed 

1=1 
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E"=i |ajf sj(^) = {^WH^H\cj)). Hence, R exists and 
is equal to R'. That is to say, Eq. (jH) is well-defined and 
its extremum is attained by picking H and \4>) so that 
M = and Aj(H)t/;ie (-7r,7r]. □ 



Appendix B: Proof of Theorem [2] and the conditions 
for equality in Eq. p6p 

Proof of Theorem [H I only need to show the second half 
of the inequality in Eq. (ITS]) as the first half follows di- 
rectly from it. More precisely, from the second half of 
Eq. m, ^{U)p < i^iUV)p+u{V-')p = ,y(UV)p+i^iV)p; 
and similarly i^{V)p < v{U-^)p + v{UV)p = v{U)p + 
v{UV)p. And from Eq. (|14ap . it suffices to prove this the- 
orem by showing its validity for all jl^^^ {j = 1, 2, . . . , n). 
Let e > be a small real number. Then, 



3,1 



(Bl) 



I now follow the strategy used in Ref. [27| to bound the 
eigenvalues of UV^ . By writing UV^ in the orthonor- 
mal basis {\4'^)}'^^i, I can regard UV^ as a matrix with 
matrix elements 



Let me consider the effect of on the non-degenerate 
eigenvalues of U first. Suppose e'^» is a non-degenerate 
eigenvalue of the unperturbed unitary matrix U . I define 

I 



an invertible matrix F by 



ifj^fc, 
= { ea if j = k = a, 

1 otherwise, 



(B3) 



where a > is a constant to be determined later. Now 
I can bound the location of the corresponding perturbed 
eigenvalue of UV^ by applying Gerschgorin circle theo- 
rem to the isospectral matrix FUV^F^^ . Specifically, 
the center Cj and radius Rj of the jth Gerschgorin circle 
on the complex plane for the matrix FUV'^F~^ are given 
by 



Cj = {FUV^F-')^iUV^)^ 



(B4) 



and 



R, = 



J2\{FUV'F-^) 



if j = a. 



(B5) 



if j ^ a. 



Combined with 



Ojk 



(B6) 



1=1 



where 5jk is the Kronecker delta and 



\ix _ 2 = 2sin — < |a;| whenever < a; < — , 

I I 2 - ' ' 2 - - 2 

(B7) 

I conclude that the distance between the centers of the 
ath and jth Gerschgorin circles with j ^ a on the com- 
plex plane obeys 



\Ca-C,\> 



> 



> 



e 3 



e-a y^^e'^"' -1 

=1 



n 

1=1 

2^|e-<-l 

11 



kU\Ay\ 



£=1 



(B8) 



And by the same argument, for any 1 < x ^ y < n, 



n 

5:(e-^r-l) (C|<^r)('/'ri0,^) 

n 

<^|e-«r _i 



(B9) 



I 

Hence, the sum of radii of the ath and jth Gerschgorin 



1=1 
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circles satisfies 



< 



- + e(n - 2) + e^a{n - 1) 
a 



E 



- + e{n - 2) + e'^an 



El 

1=1 



(BIO) 



Since e*^" is a non-degenerate eigenvalue of the matrix 
C/, I can always find a, e > satisfying 



ea < 



4n — n 



2n 



for j = 1,2 . . , n. 

Eq. (IB13|) also holds even if eigenvalues of U are de- 
generate. To show this, all I need to do is to modify the 
above proof as follows. Suppose e*^" is an r-fold degen- 
erate eigenvalue of U. And denote the r-dimensional 
Hilbert subspace making up of the corresponding degen- 
erate eigenvectors of U. For e sufficiently close to 0, is 
strictly diagonally dominant. Hence, the r x r submatrix 
y 1^ formed by retaining only the rows and columns of 
corresponding to the span of H is diagonalizable. In 
other words, there exists a basis {\(f>f)}"^i such that 
(Bll) and V'^ly^ are simultaneously diagonalized. Regarding U 
and as matrices in this basis, I know that 



and 



a > max -. ■ 

exp [i9^ ) — exp [i- 



(B12) 



With these choices of a and e, the ath Gerschgorin circle 
will be disjointed from all the other Gerschgorin circles 
because \Ca — Cj\ > Ra + Rj for all j ^ a. So, according 
to Gerschgorin circle theorem, there is exactly one eigen- 
value of UV^ located inside the circle centered at Ca and 
radius Ra < e^ck^^ |^/| = C)(e^) on the complex plane. 
Hence, if eigenvalues of U are non-degenerate, the eigen- 
values of UV^ for sufficiently small e are given by 



E' 

fe=i 



'+ee 



1 



^"1 



0(.^) 



(B13) 



•) = = Oif j^aand 

(B14) 

provided that e > is sufficiently small. I now express 
UV^ in this basis and replace the diagonal matrix F in 
Eq. jlSI by 



Fjk - 



if 
ea if j 

1 otherwise 



k and ef 



^ a J 



(B15) 



Then, if 6*^^ ^ 9^, the inequahties in Eqs. jHH and (|BT0)) 
constraining the centers and radii of the ath and jth 
Gerschgorin circles of the matrix FUV^F^^ still apply. 



Whereas in the case of 9 



u 



Ca — Cj 



n 

E(|(<^al^^'')| 

n 



By Taylor's formula with remainder. 



-2E {mm' -Mm') {mm'-m\cbi] 



\Ca-C,\= Bje + Oie^) 



where 



^E(i('^ai^^)r-i(^7i^nr)(i('^ai'^)^)r-i(<^7i^i^)r) 



1/2 



r 



(B16) 



(B17) 



(B18) 



It is important to note that Bj > if Cq 7^ Cj and that FUV^F ^ obeys 

the magnitude of the 0(6"^) remainder term in Eq. (IB17[) n n 

is less than or equal to that of the Bje term provided that Ji - < e'^a{n — r) ^ 1^/ 1 < ^^an ^ \9j\ if 9j' = 9^ . 

e is sufficiently close to 0. From Eqs. (|B14p and (jBlSp . 

the radius of the jth Gerschgorin circle of the matrix (B19) 

Suppose the set {j : Ca — Cj} has r' < r elements. Then, 
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by choosing a, e > satisfying Eq. (jBllI) . 



a > 



re I 



J : ef^ei; |exp {i0^) ~ exp {iOf) \ 



and 



4ean 



< 



j: e"=ev and B^t^O 



B 



3^ 



(B20) 



(B21) 



the r' Gerschgorin circles with a common center Ca and 
each with O(e^) radius will be disjointed from the rest 
of the Gerschgorin circles. Hence, by Gerschgorin circle 
theorem, there are exactly r' eigenvalues of UV^ each 
obeying Eq. (|B13|) . Hence, Eq. (IB13|) also holds f or the 
degenerate eigenvalue case. I also remark that Eq. (|B13[) 
resembles the Rayleigh-Schrodinger series truncated at 
the terms for time-independent perturbation of Her- 
mitian operators. 



The 



arg 



argument 

ieer \IaJJ\. 



of Eq. (|BT3)) equals e\ 



In order words. 



the arguments of the eigenvalues of UV^ obey 



9^ 



9^^ 



fe=i 



^71 



-0(6^) mod 27r. (B22) 



Since e is a small positive number and all arguments are 
written in their principle values, Eq. (jB22p implies 



3 



< 



< 



1^ 



3 I 



n 

fe=i 



O(e^) 



0(e2)(B23) 



Note that for sufficiently small e > 0, the equality in the 

|2 



first line holds if and only if 



^"1 



< 



TT. Now, applying the eigenvalue perturbation and stabil- 
ity results for the sum of two diagonalizable matrices in 
Ref . [131 , I conclude that the eigenvalues of the positive 
semi-definite Hermitian matrix 

H{U,V,e) 

n 
n 

= E [l^ll(f^) \'l^](U)){<l>jm + e\e\]{V) |0j(^))(<^j(F)| 

= Ha{U) + eHbiV) (B24) 

are precisely those given in the last line of Eq. ((B23l) . 
(Similar to the above analysis for the degenerate eigen- 
value case, this result is also true when eigenvalues of 
Ha{U) are degenerate. The key of the proof is to care- 
fully pick a basis so that Ha{U) and are simul- 
taneously diagonalizcd for every degenerate subspace of 
Ha{U).) Combined with Eq. (IH23)) . I conclude that 



m{UV')<X\{H{U,V,e)) + 0{e^) 



(B25) 



for J = 1, 2, • • • ,n. 

According to Corollary 6.6 in Ref. [l^ . 

m m 

E \0\]{UV^) < J2 ^\^^{U, V, e)) + 0{e' 



i=i 



3 = 1 



m 

<J2[xj{Ham + X]{em{V))] +0( 



m 

J2[\(^\3iU) + mj{V)]+0{e'). (B26) 



Interestingly, this inequality is the Hermitian analogy of 
what I need to prove here. 

Now, by iteratively applying Eq. (|B26p to UV = 

q terms 

, ^ 

^^1/9-^1/9... ^1/9 I get 



7n m 

Y.\(^\]{uv) <Y,[\e\j{u) + \e\j{v) 



O(-). (B27) 

q 



By taking the limit q 
equality in Eq. (HI 



-oo, I obtain the second in- 

□ 



Remark 4. Note that for each j — 1,2, ... ,n, the equal- 
ity of Eq. (B2B\) holds if and only if 0f^\ 9^ and 
are all non-negative or non-positive for all k whenever 
{4>V\4'Y) 0. And from the proof of Corollary 6.6 in 



Ref. [Is], the equality of Eq. \B20f) holds if and only if 
the spans of {\i]{H{U,V,e)))}f^^, {|4(#a(C/)))}™ i and 

{\^]{Hb{V)))}Y=i agree. Suppose the vector fl consists of 
r distinct fij 's. Then using the above two observations 
and by induction on r, one can prove the following nec- 
essary and sufficient conditions for the second inequality 
in Eq. \10fj of Theorem [H to become an equality. The 
detailed proof is left to interested readers. 

1. The n-dimensional Hilbert space T-L on which U,V 
and UV act can be written as the direct sum 
^j^iUj. Moreover U, V and UV are simul- 
taneously block diagonalizcd with respected to this 
direct sum decomposition of T-L. That is to say, 
{(l)\U\il>) — whenever and \'ip) belong to dif- 
ferent Hj 's. And similarly for V and UV . 

2. The ordering of absolute values of the arguments of 
eigenvalue of U, V and UV respects the direct sum 
decomposition of Ti in the sense that \S\^(U\^,) > 

l^lt' ( ^l-H ., ) all k, k' whenever j > j' . And sim- 
ilarly for V and UV . Furthermore, when calculat- 
ing using Eq. the absolute values of the 
arguments of the eigenvalue in each of the corre- 
sponding diagonal blocks in U, V and UV are asso- 
ciated with the same value of /ij . 

3. If the /ij corresponding to "Hj is non-zero, then 
U\-^., V\^, and UV\^, can be further simultane- 
ously block diagonalizcd with respected to the direct 
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sum decomposition "Hj = Hj ® . Furthermore, 
^t(^l«+) ^nd ei{UV\^^) > for all 
k; while el{U\^-), Ol{V\^-) and Oi{UV\^-) < 
for all k. 

Since the number of independent conditions for equality 
in Eq. il6\} increases with the dimension n of the unitary 
matrices as well as the value m used in i>{-)pim], I con- 
clude that as n or m increases, it is harder and harder 
for the equality to hold provided that the unitary matrices 
are drawn randomly from the Haar measure of U{n). I 



verify this assertion in Fig. [II which shows the plots of 
i'{UV)p[rr^] against v{U)p[^] +v{V)p[,..]- 
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